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We re-evaluate the electromagnetic corrections to rj — * 3k decays at next-to-leading order in the 
chiral expansion, arguing that effects of order e 2 (m^ — m u ) disregarded so far are not negligible 
compared to other contributions of order e 2 times a light-quark mass. Despite the appearance of 
the Coulomb pole in 77 — * n°n + n~ and cusps in 77 — * 3n°, the overall corrections remain small. 
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1. Introduction 

The decay v\ — > 3 % is particularly interesting because it is forbidden by isospin symmetry and 
thus can only happen via isospin breaking due to either strong interactions, 

J#qcd(x) = — -{dd-uu)(x) , (1.1) 
which are proportional to the light-quark mass difference — m u , or electromagnetic interactions, 

^qedW = -j J dyD^(x-y)T(j^x)j v (y)) , (1.2) 

which are proportional to the electric charge squared. A well-known low-energy theorem by 
Sutherland [jl|] states that the electromagnetic effects in r\ —> 3n are small. This decay is there- 
fore very sensitive to — m u and hence it potentially yields a particularly clean access to the 
determination of quark mass ratios. The systematic machinery that can cope with both effects 
accurately and leads to many experimentally testable predictions like Dalitz plot parameters, the 
decay width and quark mass ratios is chiral perturbation theory (ChPT) with the inclusion 

of electromagnetism In view of new and upcoming high statistics experiments [f)]-[n]] we 
reconsider the electromagnetic corrections to these observables. 

While early calculations of the strong amplitude at tree level yielded a decay width which is 
off from the experimental value by a factor of a few, Gasser and Leutwyler (GL) Jl2[ ] calculated the 
strong contributions at one-loop level and observed large unitarity collections due to strong final- 
state interactions (i.e. %% rescattering), but their result still differs from experiment by a factor of 
roughly 2. Thus strong corrections beyond one loop were studied using dispersive methods ©0], 



unitarized ChPT Q15|,|16|], and finally with a complete two-loop calculation Q17|]. All of these found 
considerable enhancement compared to the one-loop calculation. 

According to Sutherland's theorem, the electromagnetic contributions at tree level vanish. 



Baur, Kambor, and Wyler (BKW) J18Q studied corrections to Sutherland's theorem by evaluating 
the electromagnetic effects at one-loop level, but they found them to be very small. The motivation 
for reconsidering electromagnetic corrections at chiral order p 4 in the present work hinges on the 
fact that BKW neglected terms proportional to e 2 8, where 8 = m c i — m u , by arguing that these are of 
second order in isospin breaking and therefore expected to be suppressed even further. However, by 
restricting oneself to terms of the form e 2 m s and e 2 m, where m = {m u + m c /)/2, one excludes some 
of the most obvious electromagnetic effects like real- and virtual-photon contributions as well as ef- 
fects due to the charged-to-neutral pion mass difference AM 2 . = M 2 ± —M 2 = &{e 2 ), both of which 
scale as e 2 8. These mechanisms fundamentally affect the analytic structure of the amplitudes in 
question: in the charged decay channel v\ — > 7i°7t + 7i~ there is a Coulomb pole at the % + %~ thresh- 
old s = 4M 2 ± , while in the neutral decay channel Y\ — > 371° the pion mass difference induces a cusp 
behavior at this threshold due to % + %~ — > 7r°7r° rescattering, compare e.g. Ref. Jl9|]. In contrast, 
the corrections identified by BKW are all polynomials (due to counterterms) or quasi-polynomials 
(due to kaon loop effects) inside the physical region. Furthermore, Sutherland's theorem guaran- 
tees that at the soft-pion point these corrections scale as e 2 m only (and not e 2 m s ); hence the relative 
suppression of the neglected terms is of the order of5/m«2/3 and therefore not a priori small. 
The details of the presented analysis can be found in Ref. 
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2. r\ — > 3tt decay amplitudes at 0{e 2 d) 

At leading chiral order the tj — ► 37T decay amplitudes for the charged and the neutral channel, 

(n°n + n-\ri} = i(2n) 4 5 A (p Jl o + p n+ + p n - - p^^s^u) , 
(3K°\ri) = i(27i) 4 8\p n o+p 7l o+p Jz o-pr 1 )g/n(s,t,u) , 

follow from the tree graphs where the fields in the LO Lagrangian are diagonalized by use of the 
TJ7T mixing ang le £ = (V3/4) x 8/(m s — m) + ff(8 3 ). For the charged channel, the Mandel- 
stam variables s = (p^ — p n o) 2 , t = (p^ — p n + ) 2 and u = (p^ — p % - ) 2 are related by s + 1 + u = 
M 2 + M 2 + 2M 2 = 3^q. Expanding in isospin breaking parameters up to 0(8, e 2 ,e 2 8) and rewrit- 
ing everything in terms of physical observables like M 2 , M 2 , AM 2 , and the pion decay constant F n 
yields for the charged amplitude at leading chiral order 



3V3F2\ 1+ Ml-M], ]- Q 2 16V3F 2 M 2 k0 " ^ 

An additional electromagnetic term of order e 2 8 cancels the pion mass difference that is implicitly 
included in Sq. The LO charged amplitude is completely proportional to 8 and can thus be ex- 
pressed in terms of the quark mass double ratio Q 2 = (m 2 — m 2 )/ (m 2 d — m 2 ) which is particularly 
stable with respect to strong higher-order corrections [Q, |T]|. The amplitude depends linearly on 
s, and by inserting s 1 ^ it explicitly displays the Adler zero at s = 4M 2 /3. For the neutral channel, 
the Mandelstam variables s = (p^ — p n o) 2 , t = (p^ — p^o) 2 and u = {p r] — p^o) 2 are related by 
s + 1 + u = M 2 + 3M 2 = 3s'q. The neutral amplitude at leading chiral order, 

L0 _ Bo8_ 3(M 2 -M 2 )(3M 2 +M 2 ) 

V3F 2 Q 2 l6V3F 2 M 2 nQ ' ^ 

also carries an overall factor of 8, but it contains neither derivatives nor electromagnetic terms and 
hence it is just a constant. 

While for the charged decay invariance under charge conjugation implies that the amplitude 
£/ c (s,t,u) is symmetric under the exchange of t and u, the amplitude £/ n (s,t,u) for the neutral 
decay has to be symmetric under exchange of all pions and thus all Mandelstam variables due to 
Bose symmetry. All calculations of both tj — > 371 decay channels performed in ChPT so far used the 
following relation between the charged and the neutral amplitude that utilizes isospin symmetry, 

£/ n (s,t,u) = £/ c (s,t,u) -\-&Z c (t,u,s) + £/ c (u,s,t) . (2.4) 

However, this relation is only valid at leading order in isospin breaking and can not be used at 
G(e 2 8). This is most easily seen by the fact that e.g. photon loop contributions do not respect it. 



B 8 f 3(s-s c ) + 2AM 2 \ {3s — 4M 2 ) (3M 2 + M 2 , 



Furthermore, an explicit check using the LO amplitudes ( |2.2| ) and < \2.3\ ) shows that this relation is 
only valid except for terms proportional to 8 x (3sq — 3sq) = 8 x 2AM 2 . Thus both decay channels 
have to be calculated separately. 

At next-to-leading chiral order the tj — > 371 decay amplitudes receive various strong and elec- 
tromagnetic contributions. Besides renormalization effects and TJ7C° mixing at NLO, both strong 
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and electromagnetic NLO low-energy constants enter the calculation via the tree diagram with 
a vertex from the NLO Lagrangian. The 1 = rescattering of intermediate pions is particularly 
important since it gives rise to about half of the total NLO corrections, cf. Ref. [|l2|]. For the 
charged decay, the infrared (IR) divergences due to virtual-photon loops are canceled by including 
real-photon radiation (bremsstrahlung) in the soft-photon approximation at the amplitude-squared 
level. The exchange of a virtual photon in the final state between the charged pions leads to a 
triangle loop function which has some interesting features: both the real and the imaginary part 
are IR-divergent; while the IR divergence in the real part is cancelled against bremsstrahlung con- 
tributions, the imaginary part can be resummed in the Coulomb phase; furthermore it contains a 
kinematical singularity at threshold s = 4M^± , the Coulomb pole. The kinematical singularities 
in the radiative corrections, i.e. the Coulomb pole and the kinematical bremsstrahlung singularity 
at s = (M^ — Mjj.0) 2 , are part of the universal soft-photon corrections, see e.g. Ref. [22], that are 
usually already applied in the analysis of the experimental data in order to perform a meaningful 
fit of the Dalitz plot distribution. Hence we omit the unobservable Coulomb phase and subtract the 
Coulomb pole and the kinematical Bremsstrahlung singularity from the squared amplitude for the 
derivation of the numerical results. 

We have checked both the charged and the neutral amplitude in several ways: they are finite 
and renormalization-scale independent and, by use of relation (2^), both amplitudes reduce to the 
results of GL at 0(8) and of BKW [|l|] at 0(e 2 ). In addition, we have checked explicitly 
that at the soft-pion point the new electromagnetic corrections are relatively suppressed in 
comparison with the BKW corrections only by 8/m rj 2/3, as expected. 



3. Results 



It is well known that one has to go beyond one-loop chiral order to obtain a phenomenolog- 
ically successful representation of the 77 — > 37T decay amplitudes. Since we focus on the electro- 
magnetic contributions anyway, all results are given as relative corrections to the results obtained 
via the purely strong amplitude at 0(8) which corresponds precisely to the GL amplitude. Further- 
more we only consider uncertainties in the electromagnetic corrections and diregard higher-order 
hadronic corrections. Hence the purely electromagnetic corrections at 0(e 2 ) corresponding to the 
BKW amplitude as well as the new mixed corrections at 6(e 2 8) are successively added to the 
strong amplitude and their uncertainties are obtained by variation of the rather unknown electro- 
magnetic low-energy constants. 

A comparison of the different electromagnetic contributions at the amplitude level is given in 
Figs. [I] and ||, where the real and the imaginary parts of the amplitudes at different orders in isospin 
breaking (GL, BKW and DKM) are shown separately for each decay channel. The amplitudes are 
plotted along the lines with t = u, the vertical dotted lines show the limits of the corresponding 
physical regions. The BKW corrections are purely real for both decay channels, as no pion rescat- 
tering diagrams contribute at that particular order, and thus the imaginary parts of GL and BKW 
coincide with each other. Furthermore, all imaginary parts are independent of any low-energy 
constants and therefore plotted without an error range. The charged decay amplitude exhibits the 
kinematical singularities at threshold s = 4M^± (Coulomb pole & phase) that are retained here for 
illustrative reasons, while the IR divergences in the amplitude are cured by hand. By our choice 
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Figure 1: Real and imaginary parts of the charged amplitudes GL (dashed/black), BKW (dot-dashed/blue), 
and DKM (full/red) for t = u. The inserts show the region close to the two-pion thresholds. The line widths 
in the real part indicate the error bands. 
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Figure 2: Real and imaginary parts of the neutral amplitudes GL (dashed/black), BKW (dot-dashed/blue), 
and DKM (full/red) for t = u. The hatched regions denote the error bands. 



of the neutral pion mass for the isospin limit, the threshold cusp in the GL amplitude is artificially 
removed from the physical threshold energy to s = 4M^ . Since the LO neutral decay amplitude is 
constant and also at NLO the dependence on s is weak, the overall variation and thereby the scale 
of Fig. § is very small. Thus the electromagnetic uncertainties as well as the expected cusp at the 
energy of the % + %~ threshold inside the physical region at &(e 8) are clearly visible. Furthermore 
the figure shows that the size of the DKM contributions is comparable to or even larger than the 
BKW contributions. 

The Dalitz plot for the charged decay is described in terms of the symmetrized coordinates 

-T+-7L _ V3(u-t) _37p i _ 3[(M 11 -M 7[0 ) 2 -s] 

2M n Q c 



V3- 



y 



Qc 2M r] Q c ' Q c 

whereas for the neutral decay it is convenient to use the fully symmetrized coordinate 



(3.1) 



(3.2) 
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Kl 2 


a 


b 




GL 


0.0325 


-1.279 


0.396 


0.0744 


ABKW 


(-1.1±0.9)% 


(+0.6±0.1)% 


(+1.4±0.2)% 


(+1.5±0.5)% 


AD KM 


(-2.4±0.7)% 


(+0.7±0.4)% 


(+1.5±0.7)% 


(+4.4±0.4)% 



Table 1: Dalitz normalization and slopes for rj — > 7T resulting from the GL amplitude and relative 

electromagnetic corrections due to BKW and DKM amplitudes. 





Kl 2 


10 2 x a 


Z 2 /ndf 


GL 


0.269 


1.27 


= 1 


ABKW 


(-1.1±0.9 


)% (+3.7±0.5)% 


0.99 


AD KM 


(-3.3±1.8 


)% (-0.2±1.0)% 


6.20 


ADKM(w/o cusp) 


(-3.3±1.8 


)% (+5.0±1.1)% 


0.35 



Table 2: Dalitz normalization and slopes for 77 — ^ 3?r resulting from the GL amplitude and relative electro- 
magnetic corrections due to BKW and DKM amplitudes. For details, see main text. 

Here Q c = T + T+ + T- =M n - 2M K ± - M^o and Q„ = T\ + T 2 + T 3 = M n - 3M n o, where the 7} 
are the kinetic energies of the respective pions in the v\ rest frame. In order to obtain the Dalitz plot 
parameters the squared absolute values of the decay amplitudes are expanded around the center of 
the corresponding Dalitz distribution according to the standard parameterizations 

\K{x,y)\ 2 = \^Y c \ 2 {\+ay + by 2 + dx 2 + ...} , 
\<(x,y)\ 2 = \^V n \ 2 {l + 2az + ...} . 

For the charged channel, the purely strong results and the successive electromagnetic relative cor- 
rections are given in Tab. [T[ The normalization gets reduced and the slopes tend to increase. All 
corrections are at the percent level, but the BKW corrections do not represent a valid estimate of 
the dominant electromagnetic corrections. The Dalitz plot parameters for the neutral channel are 
shown in Tab. ||. By DKM(w/o cusp) we denote a fit to the inner part of the Dalitz plot with z < Zo 
such that the border region from the cusp outward is excluded, since the cusp is incompatible with 
a simple polynomial fit. As for the charged decay channel, the normalization gets reduced by a few 
percent. Here, the corrections of €?(e 2 8) are even bigger than those of 0(e 2 ). Note that the cusp 
effect leads to the single biggest modification of any Dalitz plot parameter: trying to fit the cusp 
with the polynomial function reduces a by 4% (compare ADKM to ABKW), while excluding the 
cusp region increases it again by more than 5%. The significance of this non-analytic structure is 
also reflected in the fit quality as quantified by the quoted % 2 / ndf values. This behavior can easily 
be understood qualitatively by looking at the Dalitz plot distribution at &{e 2 8) presented in Fig. |[ 
The distribution can be thought of as a bowl with three symmetrical handles where the surface 
bends down due to the cusp at constant % + %~ threshold in s, t and u each. 

The decay widths that can easily be calculated from the Dalitz plot distributions are given in 
Tab. ||. For comparison we also show the result for the charged amplitude without subtraction of the 
universal soft-photon corrections, named DKM(uc). For the neutral decay, the shifts of the width 
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Figure 3: T7 > 37T Dalitz plot distribution corresponding to the DKM amplitude. 



T] — > K°K + K t] — > 3k° 


r = T n /T c 


T GL 154.5 eV 222.8 eV 
af bkw (-i.0±0.9)% (-1.1 ±0.9)% 
af dkm (-i.9±0.5)% (-3.3 ±1.8)% 

Ar DKM(uc) (_1.0±0.5)% 


r GL 1442 

Ar BKW (-0.1 ±1.2)% 
Ar DKM (_1.4±1.8)% 
Ar DKM(uc) (_ 2 .3±1.8)% 


A gBKW (+0.24±0.22)% (+0.28 ±0.22)% 
aq vkm (+o.48±0.12)% (+0.84 ±0.46)% 
A(2 DKM(uc) ( + o.24±0.12)% 


r c /„ oc Q 4 



Table 3: Decay widths, branching ratio and quark mass ratio Q. 



are completely dominated by the corrections of the Dalitz normalizations and the DKM corrections 
are twice as large as the BKW corrections. When looking at the neutral-to-charged branching ratio, 
the difference between the BKW and the DKM result becomes even larger, since at &(e 2 ) the 
relative shifts nearly cancel in the ratio. Finally, by use of the (approximate) relation T c / n oc Q~ A 
one can read off the corrections to the quark mass ratio Q which are also quoted in Tab. ||. This 
relation does not hold for the BKW contributions, but the additional error in factorizing Q 1 2 in the 
complete amplitude can be safely neglected. In order to purify a real measurement (which includes 
electromagnetic effects) such that the purely strong value for Q can be extracted, the opposite shift 
has to be applied. Note that there are also ongoing theoretical efforts to obtain Q from a dispersive 
treatment of 77 — > 37T decays ||23|]. 



4. Summary & Outlook 

Although electromagnetic contributions to T] — > 37T decays are small in general, they ought to 
be accounted for in high-precision studies. The effects at isospin breaking 6{e 2 8) are as large as 
the effects at G[e 2 ). 



7 



Electromagnetic effects in rj — > 37Z 



Christoph Ditsche 



Since the present calculation is performed in ChPT at NLO and hence includes the rescat- 
tering effect leading to the cusp only at LO in the quark mass expansion, it is not suited to serve 
for a precision extraction of 7171 scattering lengths from the cusp effect, but rather illustrates the 
phenomenon. The theoretical framework perfectly suited for such an extraction is non-relativistic 
effective field theory [^4] -28]. The present calculation is in a sense dual to the non-relativistic 
calculation as it aims at predicting electromagnetic effects in those parts of the amplitude that are 
merely parameterized in the latter. 
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